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Abstract. Complex networks with expanding dimensions are studied, where the networks 
may be directed and weighted, and network nodes are varying in discrete time in the sense 
that some new nodes may be added and some old nodes may be removed from time to time. 
A model of such networks in computer data transmission is discussed. Each node on the 
network has fixed dimensionality, while the dimension of the whole network is defined by the 
total number of nodes. Based on the spectacular properties of data transmission on computer 
networks, some new concepts of stable and unstable networks differing from the classical 
Lyapunov stability are defined. In particular, a special unstable network model, called devil 
network, is introduced and discussed. It is further found that a variety of structures and 
connection weights affects the network stability substantially. Several criteria on stability, 
instability, and devil network are established for a rather general class of networks, where 
some conditions are actually necessary and sufficient. Mathematically, this paper makes a 
first attempt to rigorously formulate a fundamental issue of modeling discrete linear time- 
varying systems with expanding dimensions and study their basic stability property. 
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1. Introduction 

Many real-world networks appear to be different but share some similar complexity in 
such diverse aspects as varying dimensionality, intrinsic connectivity, and complicated dy- 
namics. In the endeavor of understanding the common forming mechanism of seemingly 
different networks, some recent work has already captured some essential features of various 
complex networks, particularly the small- world characteristic coined by Watts and Stro- 
gatz (1) and the scale-free degree distribution in an invariant power- law form discovered 
by Barabasi and Albert (2). It is noticeable that owing to the mathematical simplification 
for elegancy and rigor, the classical random graph theory of Erdos and Renyi (3) and the 
small-world model which inherits the same spirit have a fixed dimensionality for each given 
network. Although the scale-free model generates a growing network, the network dimen- 
sion is typically restricted to be fixed as that of a single node in order to be mathematically 
manageable when analysis comes into play. 

The main interest in this letter is to deal with a growing network with expanding di- 
mensionality, a more realistic model for the Internet to say the least. The aim is to build 
a graph framework and to lay a mathematical foundation for such a model that evolves in 
discrete time with increasing dimensions as well as complex dynamics including direction- 
ality and weights if desirable. A typical example in point is the real data transmission on 
a computer network, in which the number of nodes and the connection among them are 
both varied in time, where the connectivity may also be directed and weighted. Indeed it 
is quite interesting to think of an email network in a university or a company, where one 
computer contaminated with reinfection-enhanced virus such as the infamous W32/Sircam 
will send out many copies of the virus to computers listed in its email address book, causing 
the system server overloaded as an immediate consequence. 

A new graph model for such networks as the aforementioned computer system is es- 
tablished in this letter, where each node on the network has fixed dimensionality while the 
dimension of the whole network is defined by the total number of its nodes therefore continu- 
ously increasing. The network is characterized by a discrete linear dynamical system, where 
some new nodes will be added and some old nodes will be removed throughout the evolution- 
ary process. New concepts of stable and unstable networks are introduced, which differ from 
the classical Lyapunov stability in several aspects. In particular, a special unstable network 
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model, named devil network, is discussed. It is shown that a variety of structural and con- 
nectional properties affects the network stability substantially. Several criteria on stability, 
instability, and the devil network are finally established, actually for a rather general class 
of networks, where some conditions are necessary and sufficient. The major mathematical 
contribution of this paper is to rigorously formulate a fundamental issue of modeling dis- 
crete linear time-varying systems with expanding dimensions and study their basic stability 
theory. 

2. A graph model of simple computer networks 

Consider an isolated local-area computer network with only one server for simplicity, 
assuming that at most one PC is being added to the network at a time. In the model, 
connections among nodes are directed and the directions may vary in time, but bi-directional 
data transmissions are not permitted. 

In a real-world network, some new nodes may be added and some old nodes may be 
removed from time to time. When a node is removed from the network at some time to, one 
treats it as an isolated node starting from to- This means that those removed nodes will not 
have any connections with the other nodes in the network for all t > to, and consequently 
all the corresponding connection weights become forever. 

Assume that the computer network has rit computers, referred to as nodes, at discrete 
time i e Z+ = {Ot^o- Let Xi{t) be the difference in data amount between the input and 
output of node i, 1 < i < rit, si time t, whose absolute value \xi{t)\ is called the storage of 
node i at time t. Then 

nt 
i=l 

is the total difference between the input and output data of the whole network at time t. 
Note that this A{x){t) is also the difference between the output and input data in the server, 
referred to as the central station, at time t, and its absolute value |A(a;)(t)| is called the 
storage of the server or the whole network at time t. Only the case of finite memories is 

considered; namely, every computer and the server have a maximum allowable storage. 

Conceptually, if the amount of data stored on a computer is less than or equal to its 
maximum storage and the amount of data stored on the server does not exceed its maximum 
at some time, then the network is running well at that moment. The network is said to be 
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stable if it runs well at all times. Otherwise, if the amount of data stored on some computer 
is larger than its maximum storage at some time, whenever in the process, the network is in 
a troublesome situation since it would require the server or some other computers to share 
the extra workloads. There is another scenario that the actual storages of all computers 
are less than or equal to their maximum allowable storages, but the server is overloaded, at 
a moment. In this latter case, the server will breakdown. Both of these two cases of the 
network are referred to as being unstable. 

Mathematically, the above concepts of stabihty and instabihty are defined for the model 
as follows. Let Si be the maximum storage of node i and Mq the maximum storage of the 
server (or the whole network). 

Definition 1. A network is said to be stable if there exists a positive constant tq < Mq such 
that for aU initial point x(0) = (xi(0), 0:2(0), . . . , x„o(0))^ G R"» satisfying Ej=i \xj{0)\ < ro, 
one has \xi{t)\ < Si, I < i < rit, and |A(a;)(t)| < Mq, for all t e Z+. Otherwise, it is said 
to be unstable. In particular, the network is called a devil network if it is unstable and, 
further, for any small positive constant r < Mq there exists an initial point x{0) e R"° with 
E"=i \xj{0)\ < r such that \xi{t)\ < Si, 1 <i < n*, for all t e Z+, and \A{x){tk)\ < aMo 
for infinitely many times tk > 0, k > 1, where the constant < a < 1 is called a scaling 
parameter of the network, but |A(a;)(t^)| > Mq for infinitely many times t'^ > tk, k > 1. 

Remark 1. 

(i) The above-introduced definition of stabihty for networks is different from the classical 
Lyapunov stability for dynamical systems (4). 

(ii) In the above definition of the devil network, the constant a, < a < 1, is determined 
by some specific requirements on the network. For example, one may choose a — 1/2 
in the data transmission model of the computer network. H \xi{t)\ < Si, 1 < i < rit, 
for all t e Z+, and \A(x){tk)\ < Mq/2 for infinitely many times tk, then each computer 
runs very well at any time and the whole network works fine at all t^. But, in the 
case of |A(x)(t'^)| > Mq with > tk, the network would be in a rapidly changing 
troublesome situation (devil behaviors) after t'f, — tk, k > 1. 

(iii) The condition t'^ > tk, /c > 1, in the definition of devil network, is not restrictive. 
Since {t^j^j^ and {tfej^^i are both infinite sequences, one can easily choose suitable 
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t'f,, satisfying t'j., > tk, k' > k > 1. 

Define the dimension of the network at time t be equal to rit, the total number of 
all the nodes in the network at time t. Ignoring nonlinear factors, simply assume that 
x{t) — {xi{t),X2{t), . . . , Xnt{t))'^ e R"* satisfies the following discrete hnear system: 

x{t + l) = A{t)x{t), t e Z+, [1] 

where A{t) = (ajj(t)) is the coupling matrix of the network, which is an rif^i x rif matrix. 
Given any initial point x{0) E 11""°, the solution x{t) of system 1 can be written as 

x{t + l) = D{t)x{0), t G Z+, [2] 

where 

D{t) = A{t)A{t-l)---A{(}). [3] 

Consider the following special case of the above network model in the rest of this section, 
where the matrix A{t) has entries taken from the triple {—1, 0, 1} for all t e Z+, which is 
called a T-matrix: 

(Hi) The matrix A{t) of system 1 is an n^+i xn^ T-matrix. Assume that its entries cnjit), i ^ 
j, are evaluated at t in the following way: aij{t) = — 1 if node i sends data to node 
j; aij{t) = 1 if node i receives data from node j; and aij{t) = if there is no data 
transmission between nodes i and j. It is natural to set aii{t) — 0. Clearly, the matrix 
A{t) is antisymmetrical: 

(H2) The network initially has two nodes at t — 0, i.e., no — 2, and the number of nodes in 
the network increases by one at a time; that is, rit :— t + 2 ior t > 0. The new node 
does not send or receive any data from all the old nodes at time t; that is, 

at+3jW = 0, l<j <t + 2, teZ+. 
(H3) The matrix A{t) is in the following form: 

A{0) - ( ) ' ^(2^) = ( ) ' ^ ^ 1' 
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where 



-J B{2t + l) 
A{2t + l)^\ -B^{2t + 1) C(2t + 1) |,ieZ+, 




V 021 022 • • • 02,2t+l / 



[5] 



is a 2 X (2t + 1) T-matrix, C{2t + 1) is a (2t + 1) x (2t + 1) antisymmetrical T-matrix, 
with zero block-matrices in compatible dimensions, and 

' -1 \ 



J : = 



1 



It is clear that J is a 2 x 2 antisymmetrical matrix and satisfies = —I2, where I2 is 
the 2x2 identity matrix. 

It follows from (H3) by induction that 



D(2t)^l ^ ), D(2t+1)^ 



I 



\ 



\ 



-S^(2i+l)J 




, i e Z+. 



[6] 



For any initial point x(0) = (a, 6)^ e R^, it follows from 2 and 6 that the corresponding 
solution x{f) can be written as 



x(2t+ 1) 



( -h\ 
a 

\ ) 



; x{2t + 2) 



/ a \ 

b 

u{2t + 2) 

V y 



, t e Z+, 



[7] 



where 



u{2t + 2) = -B'^{2t + 1) Jx(0) 

= (6116 - 621a, hi2h - 622a, • • • , &i,2t+i& - 62,2t+ia)^(2t + 1). 
Since hij{2t + 1) G {-1, 0, 1}, it follows from 7 that for all i e Z+ and for all l< j < rit, 

\xj{t)\ < \a\ + \b\. 



In addition, it follows from 7 that 
A{x){2t+1) = -i) + a, 

A(x)(2t + 2) = a(l - EiLV b2^{2t + 1) 



b{l + j:tfbu{2t + l)), 



[9] 
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which imply that 

|A(x)(2t + 1)1 < \b\ + |a|, 



\A{x){2t + 2)\ < \a\l-Y:iLfb2i{2t + l) + |6| 1 + E-ll' &H(2t + 1) ,teZ+. 



2t+l 



[10] 



As noted in Remark 1, for this computer network model, one may choose the constant 
a = 1/2 in Definition 1. The discussion on the stability of system 1 with this choice is 

divided into the following two cases: 

Case I. Suppose that the two sequences { E^IV bii{2t + 1)}^°^^ and { E^IV hi{2t + 1)}^°^^^ 
are bounded. Then, for any {2t + 1) x {2t + 1) antisymmetrical T-matrix C{t), the network 
is stable. 

In fact, by the assumption there exists a positive constant 7 such that 



2t+l 

1- J2b2ii2t+1] 

i=l 



2t+l 

1+Y.bu{2t + 1) 

i=l 



< 7, teZ 



+ 



Denote 

Co := M{Sj : l<j<nt,te Z+} [11] 

and only consider the situation where Co > in the following. 
For any initial value x{0) = (a, 6)^ with 

\a\ + \b\ < min{Co, Mq, M0/7}, 

it follows from 8 and 10 that the solution x{t) satisfies 

\xjit)\<Co<Sj, l<j<nt, \A{x){t)\ < Mo, teZ+. [12] 

Hence, the network is stable. 

Case II. For any {2t+l) x {2t+l) antisymmetric T-matrix C{t), the model is a devil network 
if and only if at least one of the two sequences { EtV &ii(2t+l) 1 _^ and { Etf hii^t+l) \^_^ 
is unbounded. 

The necessity follows from the conclusion of Case I. 

To show the sufficiency, without loss of generality, suppose that | J2'i=i^ ^ij(2t + ^ is 
unbounded. Then, for any a € R with < |a| < min{Co, Mo/2}, there exist infinitely many 
tk ^ ^, k > 1, such that 



1 - Yl b2i(2tk + 1) 

i=l 



> Mo/\a\. 
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Consequently, by 8 and 9, the solution x{t) of system 1 with the initial value x{0) — (0, a)^ 
satisfies that, for alH > 0, 

\xj{t) \ = \a\ <Co< Sj, 1 < j < Tit, 
\A{x)(2t + l)\ = \a\ < Mo/2, 



and 

2tk+l 



\A{x){2tk + 2)\ = \a\ 



1 - Yl b2^{2tk + 1) 



i=l 



> Mo. 



Therefore, the model is a devil network. 

To this end, it should be noted that in the second case above, for some initial values the 
storages of the server (or the whole network) can oscillate more and more strongly as time 
evolves. 

Example. Consider the special case of B{2t + 1) = {Bi{2t + 1), 0), where is a 2 x t zero 
matrix and 

n /n. / • • • \ 

^'P( + l) = (_i _i ... _i j 

is a 2 X (t + 1) T-matrix. It is clear that B{2t + 1) satisfies the condition given in Case II 
above. So, the network is a devil network. It follows from 7 that 

x{2t + 1) = {b, -a, of, x{2t + 2) = {a,b,a, ■ ■ ■ , a, Of, t e Z+, 

where in the first relation is a {2t + 2)-dimensional zero row vector and in the second 
relation is a (t + l)-dimensional zero row vector. Hence, for all t e Z+, one has \xj{t)\ < 
max{|a|, |6|}, 1 < j < fit, and 

\A(x)(2t+l)\^\b-a\, \A(x){2t + 2)\^\b + {t + 2)a\. [13] 

Obviously, in the case of a 7^ 0, the storage of the server strongly oscillates as time 
evolves. This illustrates that the network runs quite well at some times, but will break down 
at some other times, when time is sufficiently large. 

Moreover, the matrix B{2t + 1) in this example describes the phenomenon that the 
second computer in the network continuously sends data a 7^ to each of the other t + 1 
computers. Although the burden of each of these t + 1 computers received from this second 
computer is equal to a, which is not heavy if |a| is small, it gives an extra load to the server. 
If the network is an email system, this example explains why a virus-contaminated computer 
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can cause the server to breakdown since burden is continuously building up on the server in 
this way as described by the new model. 

Remark 2. 

(1) Since the coupling matrix A{t) of the computer network discussed in this section is 

a T-matrix, all the entries of D{t) defined by 6 are integers. So, system 1 cannot 
be chaotic in the sense of Li-Yorkc. However, when the connections are weighted or 
time- varying, the linear system 1 may become chaotic in the sense of Li-Yorke (5), 
which will be further discussed elsewhere in the near future. 

(2) In the above example, ii y = A(x) is taken as an output of the system, the output 
according to 13 is chaotic in the sense of Li-Yorke. In fact, in this case, there is an 
uncountable scrambled set in the diagonal line {(a, a) : a e R} (6). 

3. Stability for a general linear model of networks 

Consider the stability of a general model of networks, i.e., its corresponding system 1 is 
linear, in which its connections may be directed and weighted, and its dimension, connectivity 
as well as weights may vary with time. 

Let rif be the number of all the nodes in the network at time t. Suppose that Xj{t) 
represents a quantity of some property V of node j at time t, 1 < j < rit, and x{t) — 
{xi{t) , X2{t) , . . . , Xnt{t))'^ e R"* satisfies the linear system 1, where A{t) — {aij{t)) is an 
rit+i X rit matrix and its entry (i) represents a weight, which is no longer restricted to the 
set of {—1, 0, 1}, with direction from node i to node j at time t. 

Similarly assume that each node i in the network has its own maximum quantity (e.g., 
storage) S'j, invariant in time, and the whole network has its own maximum quantity which 
may be infinite or varying with time, for the property V. 

The following discussion is divided into two cases: (1) the maximum quantity for prop- 
erty V of the whole network is time- invariant, which can be either finite or infinite; (2) the 
maximum quantity for property V of the whole network is time-varying. 
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3.1. Networks with time-invariant mclximum quantity of property V 



Let Mq be the maximum quantity for property V of the whole network, which is a positive 
constant or infinity. In this case, the definitions of stable, unstable, and devil networks are 
similar to those given in Definition 1 in Section 2. 

It is clear that for any given initial point x(0) e R"°, the solution x{t) of system 1 can 
also be written as 2, with D{t) — (rf^ (i))„^_^^xno satisfying 3. 

Next, the stabihty and instabihty of system 1 are studied for the case of Mq < oo. 

Theorem 1. Assume that the maximum quantities Mq and Si for property V of the whole 
network and of each node i are both finite. Then, the network described by system 1 is 
stable if and only if { Yh=i dij{t)Y_ is bounded for all 1 < j < no, and moreover there 
exists a positive constant (5 such that 



\dij{t)\ < (3Si, 1 < i < no, l<i< Ut+i, t e 



[14] 



Proof. First, the sufficiency is verified. By the assumption, there exists a constant 7 > 
such that 

nt+i 

~ <1, l<j <no, te Z+. 



nt+1 
i=l 



It follows from 2, 14, and 15 that, for all t eZ+, 



no 



no 



\xi{t + l)\<J2\dij{t)\\xj{0)\<pSiJ2\^Ml i<i<nt+u 



and 



\A{x){t + l)\ = \j^lY{j:]Lid^At)xj{0)) 

< j:t=i I E^ir di,{t)\\xjm < 7E^ii \xM\- 

So, for any initial point x(0) G R"° with J2]Li \xj{0)\ < ro, where 

ro = mm{Si, S2, . . . , Sno, 1//3, Mq, M0/7}, 



[15] 



[16] 



[17] 



one has 

\xi{t)\ <Si, l<i< nt, \A{x){t)\ <Mo,te Z+. [18] 

Hence, the network is stable. 

Then, the necessity is verified. Since the network is stable, there exists a positive constant 
ro such that for any initial point x{0) G R"" with ^2]=! \xj{0)\ < tq one has \xi{t + 1)| < 
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Si, 1 < i < n-t+i, and \A{x){t + 1)| < Mq, for all t e Z+. Given any jo, 1 < Jo < n-o, set 
Xjo(O) = To and Xj{0) — 0, 1 < j 7^ Jo < riQ. Then, it follows from 2 that, for all i e Z+, 

+ 1)1 = MijoW^jo(0)l = '^o \dijo{t)\ < Si, l<i< rit+i, 



\Aix)it + l)\^\T:iVdij,it)xjM 
which implies that, for all t e Z"*", 

\dijo{t)\ < Si/r-Q, l<i< nt+u 



To 



nt+i 



UlYdij^it) <Mo, 



dijo{t) 



i=l 



Hence, inequality 14 holds with f3 = I/tq, and | J27=i dij{t)^^ ^ is bounded for all 1 < j < 
no. The necessity is thus verified. 
Therefore, the proof is complete. 

Theorem 2. Assume that the maximum quantities Mq and Si for property of a whole 
network and each node i are both finite. Then, the network described by system 1 is a devil 
network if there exists a positive constant /3 such that 

\dij{t)\ < pSi, 1 < J < no, l<i< nt+i, t e Z+, 

and moreover there exist two time subsequences, {tk}'^=i and {t^lfcLi with — > 00 and 
t'f^ ^ 00 as k ^ 00, such that | Z]"Ii^^ o?ij(ife)}^ ^ is bounded for all 1 < j < no and 
X^j^i ^ijo(^fe)|^ ^ is unbounded for some 1 < jo < no. 

{^t +1 / \ 1 
X]j=i dij{tk)j^ ^ is bounded for all 1 < j < no, there exists a constant 

7 > such that 



Yl dij{tk) 



i=l 



< 7, /c > 1, 1 < j < no. 



YliiJi <^ijo(^fe) 1 j^. ^ is unbounded for some 1 < jo < no. Without loss of 
generality, suppose that 



X] dijo{tf. 



1=1 



00 as /c — > 00. 



Then, for any positive constant r < min{5'i, S2, ■ ■ ■ , Sn^, 1//3, Mo, q;Mo/7}, where a, < a < 
1, is the scaling parameter for system 1, there exists ko > 1 such that 



X/ dijo{tj^) 



i=l 



> Mo/r, k > ko. 
11 



Choose an initial point x{0) e R"° with Xjo(O) = r and Xj{0) — for all I < j jo < no- 
It is clear that J2]Li ki(0)| = r. With an argument similar to that used in the proof of the 
sufficiency of Theorem 1, one can easily show that the corresponding solution x{t) satisfies 

\xi{t)\ <Si, l<i<nt, te Z+, 

\A{x){tk + 1)\ < a Mo, k>l, |A(x)(t'fc + l)| > Mo, k > ko. 

Therefore, the network described by system 1 is a devil network. This completes the proof. 

In the case of Mq — oo, the following result can be easily verified by an argument similar 
to that used in the proof of Theorem 1. 

Theorem 3. Assume that the maximum quantity Si of property V for each node i in a 
network is finite and the maximum quantity for property V of the whole network is infinite. 
Then, the network described by system 1 is stable if and only if there exists a constant /? > 
such that 

\dij{t)\ < pSi, l<j<no, l<i< nt+i, t e Z+. 

The model discussed in Section 2 is revisited here based on the results obtained above. 
In this model, no — 2, and it follows from 6 that, for all t e Z+, 

\dij{t)\ < 1 < Si/Co, l<i<t + 3, J = 1, 2, 

T:iLf da{2t) = 1, Elf d,2{2t) = -1, 

E^ir dn{2t + 1) = 1 - Efil^ hi{2t + 1), Efir d^2{2t + 1) = 1 + J^Lf hu{2t + 1), 

where Cq is defined by 11. Therefore, by Theorem 1, this network is stable if and only if the 
two sequences 

{ Y.f=i bu{2t + 1)}^°^^ and {Y.f=i &2i(2t + l)}Zo arc bounded. Further, by 
Theorem 2, this network is a devil network if and only if at least one of the two sequences, 
{ Yl1=i bu{2t + 1)}~ and | E^it^ &2i(2i + 1)}°^ , is unbounded. These conclusions are the 
same as those obtained in Section 2. 

3.2. Networks with time- varying mcLximum quantity for property V 

Let M{t) be a maximum quantity of property V for the whole network at time t. 
In this case, the definitions of stable and unstable networks are also similar to those 
given in Definition 1. For convenience, they are rephrased as follows. 
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Definition 2. A network is said to be stable if there exists a positive constant ro such that 
for all initial points x{0) = (a;i(0), a;2(0), . . .,Xno)'^ e R"° with ^"=1 ki(0)l < '^o, one has 
\xi{t)\ < Si, I < i < Tit, and \A{x){t)\ < M{t) for all t E Z+. Otherwise, it is said to be 
unstable. In particular, a network is called a devil network if it is unstable and, moreover, for 
any small positive constant r there exists an initial point x{0) G R"'' with I]I=i l^i(0)l < 
such that \xi{t)\ < Si, 1 < i < ut, for all t G Z"*", \A{x){tk)\ < aM(tk) for infinitely many 
times tk, for some constant < o; < 1, and |A(a;)(fj.)| > M{t'^) for infinitely many times ij,, 
A; > 1. 

Theorem 4. Assume that the maximum quantity Si of property V for each node i in the 
network is finite, and the maximum quantity M(t) of property V for the whole network is 

finite, at any time t G Z+. Then, the network described by system 1 is stable if and only if 
there exist positive constants (3 and 7 such that 



Proof. The proof of the theorem is similar to that of Theorem 1, except the following: (i) 
in the proof of the sufficiency, 7 in 15 and 16 is replaced by ^M{t + 1), Tq in 17 is replaced 



and Mo in 18 is replaced by M{t); (ii) in the proof of the necessity, Mq is replaced by 
M{t + 1). 

This completes the proof. 

Similar to Theorem 2, the following result can be established. 

Theorem 5. Assume that the maximum quantity Si of property V for each node i in the 
network is finite, and a maximum quantity M{t) of property V for the whole network is 
finite, at any time t G Z+. Then, the network described by system 1 is a devil network if 
there exists a positive constant (3 such that 



\dij{t) \ < (3Si, 1 < J < no, l<i< rit+i, t G Z+, 
and moreover there exist two time subsequences, {tfc}^^ and {t'i,}'^i with tk ^ 00 and 



dij{t)\ <f3Si, l<i< nt+i, Yl dijit) < jM{t + 1), 1 < J < no, t G Z+. 



i=l 



by 



ro = mm{Si, S2, . . . , M(0), 1/7}, 



— > 00 as A; — > 00, such that 



and {M-i(t', + l)Ei\"'rf..o(4)} 




is bounded for all 1 < j < no 
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Proof. The proof of the theorem is similar to that of Theorem 2. So the details are omitted. 
Remark 3. 

(1) If the number of nodes in a network does not vary; that is, the dimension rit of the 
corresponding system 1 is time- invariant, then system 1 is a classical time- varying 
discrete linear system, which surprisingly can be chaotic in the sense of Li-Yorke (5). 

(2) If the effects of internal and external nonlinearities on a network are considered, the 
corresponding system is by nature nonlinear, for which the stability and complex dy- 
namical behaviors need to be further addressed in the future. 

Remark 4. In the present letter, only the simplest possible model of a localized isolated 
computer network with one server is considered. A realistic computer model, however, has 
more than one server in general, which becomes more mathematically involved, leaving a 
challenging topic for future research. 
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